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^ i Abstract 

New observable integer-valued numbers of the topological origin 
. were revealed by the present authors studying the conductivity theory 

j. I of single crystal 3D normal metals in the reasonably strong magnetic 

' field {B < W^Tl). Our investigation is based on the study of dynami- 

■ cal systems on Fermi surfaces for the motion of semi-classical electron 
in magnetic field. All possible asymptotic regimes are also found for 

' i3 — > cxD based on the topological classification of trajectories. 

■ 1 Introduction. 

o 

^ ' This is the survey article dedicated to the topological approach to the con- 

ductivity phenomena in metals in the presence of the rather strong magnetic 
X ■ field B. As can be shown, in the case of rather strong magnetic fields the 

■ topology of the Fermi surface plays the main role in the asymptotic behav- 
ior of conductivity a'^^ when B ^ oo. In this article we represent mainly 
the results in this area connected with the "Topological Quantum charac- 
teristics" introduced recently by the present authors and the classification 
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results for the regimes of conductivity behavior for B ^ oo based on the 
topological results in the theory of special dynamical systems. The main 
consideration is made here for the case of the "generic" (or " Topologically 
regular") open electron trajectories according to classification recently intro- 
duced in the works of the present authors (see [201 EI])- As was shown 
(|2Uj). this type of trajectories always leads to observable topological charac- 
teristics of the Fermi surface having the form of "stability zones" on the unit 
sphere (parameterizing directions of B) and the triples of the integer num- 
bers observable in every "stability zone" in the conductivity measurements. 
Let us say that this type of trajectories corresponds to the "generic case" 
and it is the only type of non-closed trajectories stable with respect to the 
small perturbations. Besides that, the trajectories of this type appear with 
probability 1 in all cases when the non-closed quasiclassical electron can be 
observed on the fixed Fermi surface. 

However, it can be shown that also other types of open orbits correspond- 
ing to more complicated "chaotic" behavior can be observed on the compli- 
cated Fermi surfaces. In this case, the direction of magnetic field should be 
chosen specially in the experiment and the behavior of conductivity reveals 
much more complicated features. Let us mention also, that the geometric 
forms of electron trajectories for such systems can be obtained as the in- 
tersections of the Fermi surface by the planes orthogonal to the magnetic 
field and in some sense these systems are all "analytically integrable" on the 
universal covering M^. However, after the identification of equivalent vectors 
modulo the reciprocal lattice these systems may become topologically very 
complicated since the two points in different by the reciprocal lattice vec- 
tor represent actually the same physical state of electron. For the generic 
irrational directions of B the global geometry of these intersections is apriori 
unpredictable and can be rather complicated as we will see below. 

Let us give here the brief historical survey concerning the "geometric 
strong magnetic field limit" in metals. 

The investigation of the geometric effects in the magnetoresistance behav- 
ior in normal metals was started first in the school of I.M.Lifshitz (I.M.Lifshitz, 
M.Ya.Azbel, M.I.Kaganov, V.G.Peschansky) in 1950's. Such, in the first pa- 
per IJ the crucial difference in the asymptotic behavior [B — »• oo) of conduc- 
tivity in the cases of closed (compact) and open periodic electron trajectories 
on the Fermi surface was pointed out. The strong anisotropy of conductiv- 
ity in the last case makes possible the observation of this phenomenon and 
the mean direction of the periodic open orbits in the experiment. The dif- 
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ferent examples of the complicated Fermi surfaces and the corresponding 
non-closed electron orbits were considered in the works [21 El- In particular, 
the open orbits, more general then just periodic ones were discovered in [2]. 
The geometry of such trajectories still was not very complicated and the 
corresponding electron motion was still "in average" along the straight line 
in the p-space. Again the strong anisotropy of the conductivity tensor was 
expected in this case and the form of a^'' used in [21 E] was actually the same 
as for the periodic open orbits. In the paper [Sj the important set of ana- 
lytical dispersion relations was considered and the existence of open orbits 
in the different parts of set parameters was discovered. Let us give here also 
the references on papers jH El El El (see also the survey articles jSl El and the 
book [TU]) where the different question connected with the geometry of open 
orbits for concrete Fermi surfaces of real metals were considered. 

The general problem of classification of different trajectories arising on 
arbitrary periodic smooth surface in 3D space as the intersections with ar- 
bitrary planes in (S.P.Novikov problem) was set by S.P.Novikov in J21 
and then was considered in his topological school at the Moscow State Uni- 
versity (A.V.Zorich, I.A.Dynnikov, S.P.Tsarev, 1980-2000). ^ Let us men- 
tion here that the most important breakthroughs were made in the papers 

IT^ where the important theorems about the non-compact trajectories 
were proved. Using this results the concept of the "Topological quantum 
numbers" observable in the conductivity and characterizing the generic "non- 
trivial" conductivity behavior in metals was introduced by the present au- 
thors in [201 ■ I^ particular, it was shown that these characteristics arise 
always when the stable (with respect to the small variations of B) behav- 
ior of conductivity different from the "simple" behavior corresponding to 
just compact electron orbits on the Fermi surface is observed. The corre- 
sponding "stability zones" on the unit sphere and the "Triples of integer 
numbers" (" Topological Quantum characteristics" ) give the non-trivial geo- 
metric characteristics of the complicated Fermi surfaces in metals. However, 
as we already said above the conductivity can reveal also some "degener- 
ate" (unstable) behavior in the special cases of rather exotic electron tra- 

"'^The general problem of S.P. Novikov has actually more general form and is connected 
with the global geometry of level curves of quasiperiodic functions with n quasiperiods on 
the plane. Let us say, however, that the case n > 3 was started to investigate recently (see 
|29j) and is still not so well studied as the case n — 3. Let us make also the reference on the 
work 'i52| where the connection of general Novikov problem with the specially modulated 
2D electron gas was considered. 
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jectories on the Fermi surface. The existence of such trajectories as well as 
the corresponding conductivity behavior was studied later in several works 
( [211122123 1211 )• Let say here that the full classification of all types of orbits 
for arbitrary smooth periodic surfaces in M.^ is finished in general now (see 
[221I2Z1) and the total picture of different asymptotic behaviors of conductiv- 
ity (as well as the conditions for each asymptotic regime) can be described 
in main order of 1/B in general case (j2Sl E])- Let us say also that the 
general topological investigation required rather large set of methods of both 
3-dimensional topology and the dynamical systems theory and we are going 
just to formulate here the corresponding statements needed for our purposes. 



2 Some remarks on the Kinetic theory. 

Let us describe briefly the standard approach to the conductivity in crys- 
tals based on one-particle quantum mechanical consideration. ^ We should 
consider the Shrodinger equation 

Atp + V(x)V^ = 

2m 

in the periodic potential l^(x) formed by the crystal lattice L: 

y(x + mill + m2\2 + m-sls) = ^(x) 

for all (mi,m2,m3) e 1?. 

The potential ^(x) has also all the rotational symmetries of the crystal 
and in all physical cases we have the symmetry x — x after the appropriate 
choice of the origin in x-space. 

The reciprocal lattice F* is generated by the vectors: 

. I2 X I3 fc I3 X ll fc ^ 

Ki = n- — : — — , So = n , = n- 



(Ii, 12,13) ' (ll, 12,13) ' " (ll, 12,13) 

in the space of electron momenta. 

The corresponding solutions of the stationary Shrodinger equation can be 
represented as Bloch functions 

^Let us note that we use here the Kinetic theory on the physical level of rigorousness. 
At the same time all use of differential topology is based on the mathematically rigorous 
topological results. 
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V'p,s(x) = e </'p,s(x) 

where <^p^s(x) has the same periodicity as the potential ^(x). 

The electron states can be parameterized then by the discrete index s 
(energy band) and a continuous parameter p = {pi,P2,P3) such that the 
states parameterized by p and p' (at the same s) different by any reciprocal 
lattice vector: 

nigi + n2g2 + n3g3 , (^1,^2,713) e 

are physically the same. 

The electron energy e at fixed s becomes now the function of p, e = es{p), 
periodic in p-space with periods equal to the reciprocal lattice vectors: 

es{p + nigi + n2g2 + n3g3) = es(p) , {rii, n2, ns) e 

The full set of functions es(p) inherits also the rotational symmetry of 
the potential V{'x) and for any function €s{p) we have a symmetry 

es(p) = es(-p) 

after the appropriate choice of the initial point in p-space. 

From topological point of view the correct phase space for any fixed energy 
band is the three-dimensional compact torus (Brillouin zone) = M.^/G 
rather then the open three-dimensional space M.^. The functions es{p) become 
then the (smooth) functions on with the same property es(p) = ^s{—p) 
for the appropriate initial point in T^. 

The collective electron structure can be described by the one-particle 
distribution functions fs{p) satisfying to Boltzmann equation. 

According to the classical results in the absence of external fields the 
functions /s(p) are given by the Fermi distribution 

fsip) = ^uqT^ (1) 

1 + e T 

for the fixed temperature T. 

Let us mention also that every function es(p) is bounded by some minimal 
and maximal values e^*", e^"^ being the smooth function on the compact 
manifold. 
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For the case of normal metals we will always have a situation T <^ ep and 
the distribution function will change only in narrow energy interval (~ T) 
near the Fermi energy es(p) ~ ei?, being equal to 1 for — es(p) 3> T and to 
for es(p) — eF 3> T. So, the energy bands with e^"^ < cf will be completely 
filled in this case (/s(p) = 1) and the energy bands with e'™" > ep will be 
completely empty (/s(p) = 0). 

Geometrically we have occupied all the states "inside" the surface es(p) = 

in 3D torus and empty all the states "outside". The union of surfaces 
es(p) = ep for all s such that e™" < < e™""^ gives then the total Fermi 
surface Sp of the given metal. 

The Fermi surface will be called non-singular if Ves(p) 7^ for all s on 
the Fermi level e^, otherwise we will call it singular. 

It will be convenient to consider also the natural covering Sp over the 
Fermi surface in which is a smooth 3-periodic surface (for non-singular 
Sp) in the p-space with the periods equal to the reciprocal lattice vectors. 

The gradient of function es{^) is the group velocity of a particle in the 
corresponding quantum state and describes the mean velocity of the localized 
wave packet formed by the wave functions '?/'s,p'(^) with values p' close to 
p. Easy to see that the group velocity is the odd function of p (after the 
appropriate choice of the initial point) 

Vs5r(-P) = -V, gr{V>) 

since all the functions es(p) are even. The total flux density 



i = 2 / . . . / V. ,.(p) /.(p)^ = 2 ^ / . . . / vmp) /.(p)^ 

as well as the electric currents density j = ei are zero for the equilibrium 
functions given by (^. In our further consideration we will deal with the 
small deviations from the Fermi distributions caused by small external 
forces (electric field) and the correspondent small response in the total electric 
current density. 

Let us put now Vgf.(p) = ^sgrijp) oil the corresponding components of 
the Fermi surface and consider just one function Vgr(p) on 5*^. 

We will use now the fact that all the components of Sp give the in- 
dependent contributions to conductivity tensor o"*'^. This means that we 
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can actually consider only one conductivity zone with the dispersion rela- 
tion e(p) = es(p) and the corresponding distribution function /(p) = /s(p) 
defined everywhere in T^. The corresponding contributions to the electric 
conductivity (in the presence of magnetic field) should be then just added in 
the 3-dimensional tensor cr*'^. 

So we will consider now only one dispersion relation e(p) and the function 
/(p) which satisfies to the Boltzmann equation for the spatially homogeneous 
case (we omit here the spin dependence of the distribution function for sim- 
plicity): 

f + F^^*|^ = /[/](p) (2) 

Here F^"^* is the homogeneous external force and /[/](p) is the collision 
integral. For our situation of low temperatures (T ~ IK) only the scattering 
on the impurities will play the main role so we will write here the functional 
/[/](p) in the general form: 



/[/](P) = /•••/ [/(P')(1-/(P))-/(P)(1-/(P'))]^(P,P') 

[/(p')-/(p)]a(p,p' 



(27rn)3 

\2TihY 

Using of the Boltzmann equation implies the applicability of the semi- 
classical approach for the electron motion and the evolution of the quantum 
state within the Brillouin zone can be described just by the "classical" dy- 
namical system: 

p = -p^'-^i 

with the only difference that p now is a quasimomentum and belongs to 
3-dimensional torus instead of the Euclidean space. 
In our case we will have 

^ext ^ ^ [^^^ X B] + eE = [Ve(p) x B] + eE 

where B is the homogeneous magnetic field and E is the small electric field 
responsible for the electric current. 
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The applicability of semiclassical approach requires then that the fields 
B and E are small with respect to the internal fields in the crystal. We will 
require also the condition -C ep where lob = eB/m*c is the cyclotron 
frequency which will permit to consider the quantization of levels as a small 
effect with respect to the classical motion picture. The electric field is going 
to be infinitesimally small so we don't put any conditions on it. Let us just 
mention here that these semiclassical conditions are satisfied very well for 
normal metals in all experimentally available fields B (the theoretical limit 
is S ~ 10^ — W^Tl). Thus we use the Boltzmann equation Q to describe 
the main part of the conductivity tensor o"''^ while the quantum phenomena 
will be considered as the small corrections in our situation. 

Now the only reminiscent of the quantization will be the dispersion rela- 
tion e(p) and the corresponding dynamical system on given by 

p = - [Ve(p) X B] + eE 

c 

For the stationary distribution function /(p) we can write now the Boltz- 
mann equation as 

^ [V6(p) X B] V/(p) + eEV/(p) = /[/](p) (3) 

Any Fermi distribution function (P) satisfies the condition /[/o](p) = 
and gives an equilibrium distribution for given temperature T in the absence 
of external fields. 

Let us now consider in details the main dynamical system: 

p = ^ [Ve(p) X B] (4) 

on for general periodic e(p). 

System Q is Hamiltonian with respect to the (non-standard) Poisson 
bracket 

6 6 6 

{Pl,P2} = -B3 , {p2,P3} = -Bi , {p3,Pl} = -B2 (5) 

c c c 

with the Hamiltonian function e(p). The bracket (0) is degenerate with the 
Casimir 

3 

C = J2P^B^ = ip-B) 
1=1 
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However, the Casimir function C is a multi-valued function on the three- 
dimensional torus and should be actually considered as a 1-form in it. 
Geometrically, the trajectories of (j3)) in will be given on every energy 
level e(p) = const by the level curves of this 1-form restricted on these two- 
dimensional surfaces. So we have locally the analytic integrability of the 
system (jH). However, the global geometry of the trajectories of (jD) can be 
highly nontrivial because of the non-uniqueness of values of Casimir function 
C. 

The picture becomes more visible if we consider the corresponding cover- 
ing defined by the reciprocal lattice in p-space. The corresponding 
function e(p) is then the three-periodic function in M?. The Fermi surface 
(the covering over the Fermi surface in T'^) also becomes the three-periodic 
(smooth) surface in the p-space. The Casimir function C is now well defined 
function in M? (height function) and the trajectories of Q on the Fermi level 
will be given by the intersections of the three-periodic surface e(p) = ep with 
all different planes orthogonal to B. The most important thing for us will 
be the global geometry of the trajectories in the planes H(B) orthogonal to 
the magnetic field. 

The dynamical system Q conserves also the volume element d?p in 
and does not change at all the Fermi distribution ((T)). So, in the absence of 
the electric field E we will have the electron distribution unchanged (up to 
the quantum corrections). Nevertheless, the response of this system to small 
perturbations will be completely different from the case B = and depend 
strongly on the geometry of trajectories of the dynamical system (0)). 

Let us now come back to the Boltzmann equation Q and analyze the 
effects of infinitesimally small electric field E. 

We will use for simplicity the r-approximation for the collision integral 
throughout this paper and just put for small perturbations 5f: 

I [h + 5f] (p) - / [h] (p) = Jt (Sf) (p) = --5/(p) 

T 

Let us make now also some (standard) remarks about r for our case of 
normal metals at low temperatures (~ IK). As it is well-known for the 
low temperatures both the electron-phonon and electron-electron scattering 
decrease (as and respectively) and the only reminiscent at T — > is 
the scattering by the impurities which is constant for T ^ 0. For the elastic 
scattering by impurities we can usually write the collision integral in the form 
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I if] (P) = f I . . . I |^pp,r [/(P') - /(P)] S{e{p) - e{p'^ ^'^ 



{27rny 



where rii is the concentration of impurities and Vpp' is the amphtude of scat- 
tering p — > p' on the impurity. Such we can put 

^[/](p) = ^n. / / kpp'r[/(p')-/(p)] 



^ J Jsp V(P') 



on the Fermi surface. 

Easy to see then that we have the property /[/] = for any function 
/ depending only on energy: /(p) = /(e(p)). So for this form of /[/] the 
corresponding operator Jt will have a lot of zero modes corresponding to the 
perturbations 5f{p) — 5/(e(p)). Practically this means that the processes of 
electron-electron and electron-phonon scattering responsible for the mixing 
of different energy levels become very small at low temperatures. The corre- 
sponding "energy relaxation time" T ^ and the characteristic 
times of energy relaxation become very big at low T. 

However the energy relaxation time Ten will not be important for us in 
the consideration of the conductivity phenomena. Indeed, the mean electric 
flux density given by 

Vgr(p)5/(p)- 



(27r/l)3 

remains zero for any distribution perturbation of form Sf{p) — 5/(e(p)) 
because of the property Vgr{—p) — — Vgj.(p). We have so that all the zero 
modes of do not give any contribution to the conductivity and we can 
omit them in our consideration. 

We will consider now only the perturbations 5/(p) such that the total 
number of particles on each energy level remains unchanged w.r.t. the initial 
Fermi distribution. So we put an additional restriction 

5/(p)-^^0 (6) 

e{p)=const ^gr{P) 

for all the energy levels. Let us say also that the only factor which can disturb 
this condition according to the full Boltzmann equation 
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is the electric field E. However, because of the equality 

dfx d^p f [ , .dfr d^p _ 

this effect has the order E"^ and should be omitted in the computing a linear 
approximation for the infinitesimally small E. We can now forget about the 
zero modes of the functional Jt using the functional space defined by ^ 
and put r to be the characteristic relaxation time on each energy level due to 
the scattering on the impurities. This time r = Tmom is then the momentum 
relaxation time and it is this quantity which is responsible for conductivity 
at T ^ 0. We just use then the approximation 

Jt [6f] (p) = --6f{p) 

T 

on the space of function satisfying (jH)) where r = Tmom has the meaning of the 
free electron motion time before the scattering by the impurity and remains 
constant at T ^ 0. 

The Boltzmann equation Q can be rewritten now in linear order of E as 

^ [Ve(p) X B] V/i(p) + eEv,,(p)^ = -^A(p) (7) 

where /i(p) is a linear in E perturbation of Jt- 

Let us introduce the vector field ^(p) on according to our system (jl)): 

e(p) = ^ [Ve(p) X B] 
We can then rewrite the equation ((Tj) on /i in the form: 

V€/i(p) + ^/i(p) = -eEv,,(p)^ (8) 

The equation (jH)) can be now solved separately on every trajectory of 
the system Indeed, let us introduce the coordinate t (time) along the 
trajectories of (^. We have then 
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on every trajectory. We put here /i = fi(t), v^^ = ^gr{t) and the value 
dfr/de is constant on the trajectories of (jD). 

The solution of (jH)) can be written now as ([TU]): 



J e^Evgrit - t')dt' (9) 



dt 

Let us introduce the "geometric" time s = teB/c where t is the time 
along the trajectories according to the system (j3]). For the derivative of p 
with respect to s we then will have the system: 

^ = [V6(p) X n] (10) 

where n = B/S is a unit vector along B. Let us introduce the notation 
Ygr{p,t) = Ygr{p(t)) whcrc p(t) IS & solutiou of the system (jH) with the 
initial data p(0) = p. The formula © can be then written as 

/i(p) = e (-^) Ev,,(p, -t)e-'/-dt (11) 

The parameter s plays the role of geometric parameter along the trajec- 
tories, such that dl/ds = v^^(s) where dl is the length element in the p-space 
and Vg^{s) is the length of component of Vgr(s) orthogonal to B. The size of 
the Brillouin zone corresponds then to s ^ pp/vp ^ m* where m* is some 
parameter with mass dimension playing the role of "typical" electron mass 
on the Fermi surface for a given dispersion relation. 

We will be interested in the " strong geometric limit" when the parameter 
m*c/eBT is much less then 1. This condition can be also written as ujbt S> 1 
where ub = eB/m*c is the cyclotron frequency corresponding to the mass m*. 
It is easy to see from (lllj) that this case corresponds to "long integrations" of 
Vgr(p) along the trajectories of (jlUj) (much longer then one Brillouin zone for 
open trajectories) such that the global geometry of the trajectories becomes 
the most important factor for the linear response /i(p). 

The correspondent current density j has the form 

Vgr(p)/l(p)- ^ ^ 
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dfr 
de 



- r+oo -| j3 

/ Ev,,(p,-t)e-*/-rft — 
Jo J (^^ 



and we have for the conductivity tensor a ([lOj, §27): 



T 



de 



<(P) 



p, -t)e-''^dt 



Easy to prove also the Onsager relations for tensor (T*'^(B): 



[12) 



It is easy to see that the expression 

■f)«;»<(P,-<)(^ 

is invariant under the action of the dynamical system We have so: 



■^)<(p)<(p,-t) 



-"^dt 



^^^^ vUpMrip) 



de 



{2Tihf 



'''^dt 



(13) 



Using now the fact that f*^(p,t) coincides with the value f^r(P) ~t) de- 
fined for the system (j3]) for — B we get the Onsager identity. 



3 Topological considerations and theorems. 

Before starting the investigation of the dependence of a^^ on the geometry 
of orbits of Q let us make here some basic definitions and formulate the 
topological results which we are going to use in our considerations. 

First of all, we consider the intersections of the smooth 3-periodic Fermi 
surface in p-space, with the set of parallel planes orthogonal to magnetic field 
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Figure 1: The level curves of the function e(p)|n in the plane close to the 
local minimum, the saddle-point and the local maximum of e(p)|n- 

B. We can consider the corresponding trajectories also as the level curves 
of the restrictions e(p)|n of the dispersion relation to the planes orthogonal 
to B which will then be the quasiperiodic functions in this planes with 3 
quasiperiods according to the standard definition. Let us make now the 
generic assumption that the restrictions e(p)|n are the Morse functions in 
any plane 11 orthogonal to B, i.e. all the critical points of this functions 
can be the non-degenerate local minima, non-degenerate saddle points or the 
non-degenerate local maxima in the plane. On the Fig. ^we show the local 
behavior of the trajectories close to these critical points and we don't admit 
here any other degeneration of trajectories in the plane. 
Let us give here the following definitions: 

Definition 1. We call the trajectory non-singular if it is not adjacent 
to the critical (saddle) point of the function V{r). The trajectories adjacent 
to the critical points as well as the critical points themselves we call singular 
trajectories (see Fig. Q^. 

Definition 2. We call the non-singular trajectory compact if it is closed 
on the plane. We call the non-singular trajectory open if it is unbounded in 
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Figure 2: The singular, compact and open non-singular quasiclassical tra- 
jectories. The signs " + " and " — " show the regions of larger and smaller 
values of e(p)|n respectively. 

The examples of singular, compact and open non-singular trajectories are 
shown on the Fig. |2l a-c. 

It is easy to see also that the singular trajectories have the measure zero 
among all the trajectories on the plane. 

In our consideration we will not make any difference between the com- 
pact trajectories of different forms and consider the global geometry of the 
non-compact trajectories only. According to the topological results on the 
classification of non-compact orbits it will be convenient to use the following 
definition concerning the global behavior of open trajectories in the planes 
orthogonal to B: 

Definition 3. We call the open trajectory topologically regular (corre- 
sponding to "topologically integrable" case) if it lies within the straight line 
of finite width in and passes through it from —oo to oo (see Fig. a). 
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Figure 3: " Topologically regular" (a) and "chaotic" (b) open trajectories in 
the plane 11 orthogonal to B. 

All other open trajectories we will call chaotic (Fig. 6j. 

We will see below that this approach to the classification of electron 
trajectories is closely connected with the form of the Fermi surface itself 
and the forms of the "carriers of non-compact trajectories" belonging to the 
Fermi surface. Thus the "topological integrability" is actually the property 
characterizing the topology of the carriers of open trajectories in (and M^) 
and has the direct analog in the theory of integrable systems although the 
reasons of " integrability" are completely different in this situation and have 
completely topological origin. 

Let us give now the basic definitions concerning the topology of the Fermi 
surface and later the topology of "carriers of open trajectories" on the Fermi 
surface for the fixed direction of B. 

Definition 4. 

I) Genus. 

Let us now come back to the original phase space = IR^/F*. Every com- 
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a) b) c) 

Figure 4: The abstract surfaces with genuses 0, 1 and 2 respectively. 

ponent of the Fermi surface becomes then the smooth orientable 2- dimensional 
surface embedded in T^. We can then introduce the standard genus of every 
component of the Fermi surface g = 0, 1, 2, ... according to standard topologi- 
cal classification depending on if this component is topological sphere, torus, 
sphere with two holes, etc ... (see Fig. |^. 
II) Topological Rank. 

Let us introduce the Topological Rank r as the characteristic of the em- 
bedding of the Fermi surface in T^. It's much more convenient in this case 
to come back to the total p-space and consider the connected components of 
the three-periodic surface in M^. 

1 ) The Fermi surface has Rank if every its connected component can be 
bounded by a sphere of finite radius. 

2) The Fermi surface has Rank 1 if every its connected component can 
be bounded by the periodic cylinder of finite radius and there are components 
which can not be bounded by the sphere. 

3) The Fermi surface has Rank 2 if every its connected component can 
be bounded by two parallel (integral) planes in M? and there are components 
which can not be bounded by cylinder. 

4 ) The Fermi surface has Rank 3 if it contains components which can not 
be bounded by two parallel planes in M^. 

On the Fig. |3]we show the topological form of the surfaces of genuses 0, 
1 and 2 (without the embedding in T^). It is easy to generalize this picture 
to the arbitrary genus g. 

The pictures on Fig. a-d represent the pieces of the Fermi surfaces in 

with the Topological Ranks 0, 1, 2 and 3 respectively. 
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It is easy to see also that the topological Rank coincides with the max- 
imal Rank of the image of mapping 7ri(S'*) vri(T^) for all the connected 
components of the Fermi surface. 

As can be seen the genuses of the surfaces represented on the Fig. 
a-d are also equal to 0, 1, 2 and 3 respectively. However, the genus and 
the Topological Rank are not necessary equal to each other in the general 
situation. 

Let us discuss briefly the connection between the genus and the Topo- 
logical Rank since this will play the crucial role in further consideration. It 
is easy to see that the Topological Rank of the sphere can be only zero and 
the Fermi surface consists in this case of the infinite set of the periodically 
repeated spheres §^ in M^. 

The Topological Rank of the torus can take three values r = 0, r = 1 
and r = 2. 

Indeed, it is easy to see that all the three cases of periodically repeated 
tori in R^, periodically repeated "warped" integral cylinders and the peri- 
odically repeated "warped" integral planes give the topological 2-dimensional 
tori in after the factorization (see Fig. E)). 

Let us note here that we call the cylinder in integral if it's axis is 
parallel to some vector of the reciprocal lattice, while the plane in is 
called integral if it is generated by some two reciprocal lattice vectors. The 
case r = 2, however, has an important difference from the cases r = and 
r = 1. The matter is that the plane in is not homological to zero in 
(i.e. does not restrict any domain of "lower energies") after the factorization. 
We can conclude so that if these planes appear as the connected components 
of the physical Fermi surface they should always come in pairs, 11+ and n_, 
which are parallel to each other in M^. The factorization of 11+ and n_ gives 
then the two tori T^, with the opposite homological classes in after 
the factorization. The space between the 11+ and n_ in can now be taken 
as the domain of lower (or higher) energies and the disjoint union 11+ U n_ 
will correspond to the union U homological to zero in T^. 

It can be shown that the Topological Rank of any component of genus 2 
can not exceed 2 also. The example of the corresponding immersion of such 
component with maximal Rank is shown at Fig. c and represents the two 
parallel planes connected by cylinders. 

At last we say that the Topological Rank of the components with genus 
g > 3 can take any value r = 0, 1, 2, 3. 

Let us give also the definitions of "rationality" and "irrationality" of the 
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Figure 6: The periodically repeated tori T^, periodically repeated "warped" 
integral cylinders and the periodically repeated "warped" integral planes in 
M3. 
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direction of B. 



Definition 5. 

Let {gi,g2,g3} be the basis of the reciprocal lattice T* . Then: 

1 ) The direction of B is rational ( or has irrationality 1) if the numbers 
(B,gi), (B,g2), (B,g3) are proportional to each other with rational coeffi- 
cients. 

2) The direction ofB has irrationality 2 if the numbers (B,gi), (B,g2), 
(B, ga) generate the linear space of dimension 2 over Q. 

3) The direction o/B has irrationality 3 if the numbers (B,gi), (B,g2), 
(B, ga) are linearly independent over Q. 

The conditions (l)-(3) can be formulated also as if the plane n(B) or- 
thogonal to B contains two linearly independent reciprocal lattice vectors, 
just one linearly independent reciprocal lattice vector or no reciprocal lattice 
vectors at all respectively. 

It can be seen also that if {li, 12,13} is the basis of the original lattice 
in x-space then the irrationality of the direction of B will be given by the 
dimension of the vector space generated by numbers 

(B, 12,13) , (B,l3,li) , (B,li,l2) 

over Q. Easy to see that these numbers have the meanings of the the magnetic 
fluxes through the faces of elementary lattice cell. 

Let us discuss now the connection between the geometry of the non- 
singular electron orbits and the topological properties of the Fermi surface. 
We will briefly consider here the simple cases of Fermi surfaces of Rank 0, 1 
and 2 and come then to our basic case of general Fermi surfaces having the 
maximal rank r = 3. We have then the following situations: 

1) The Fermi surface has Topological Rank 0. 

Easy to see that in this simplest case all the components of the Fermi 
surface are compact (Fig. [71 a) in and there is no open trajectories at all. 

2) The Fermi surface has Topological Rank 1. 

In this case we can have both open and compact electron trajectories. 
However the open trajectories (if they exist) should be quite simple in this 
case. They can arise only if the magnetic field is orthogonal to the mean 
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Figure 7: The simplest compact and open periodic trajectories. 



direction of one of the components of Rank 1 (periodic cyhnder) and are 
periodic with the same integer mean direction (Fig. [7J b). There is only the 
finite number of possible mean directions of open orbits in this case and a 
finite "net" of one-dimensional curves on the unit sphere giving the directions 
of B corresponding to the open orbits. In some special points we can have the 
trajectories with different mean directions lying in different parallel planes 
orthogonal to B. Easy to see that in this case the direction of B should be 
purely rational such that the orthogonal plane n(B) contains two different 
reciprocal lattice vectors. It is evident also that there is only the finite 
number of such directions of B clearly determined by the mean directions of 
the components of Rank 1. Let us mention also that the existence of open 
orbits is not necessary here even for B orthogonal to the mean direction of 
some component of Rank 1 as can be seen from the example of the "helix" 
represented on Fig. |Hl a. Easy to see also the this type of trajectories 
corresponds to topologically integrable case according to the Definition 3 
and gives the simplest example of topologically regular open orbit in the 
plane n(B). 

Let us point out here that both cases of compact and periodic open tra- 
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Figure 8: (a) Connected component of Rank 1 having the form of "hehx". 
Open orbits are absent for any direction of B. (b) The example of the 
Fermi surface of Rank 2 containing two components with different integral 
directions. 
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jectories were considered in details in the work pQ where the form of the 
conductivity tensor in the hmit B ^ oo was obtained for these situations. 
Let us say also, that these two situations are actually the only two cases 
where the "dynamical" coordinates {e,pz,t) considered in P (corresponding 
to the dynamical system can be introduced globally on the Fermi sur- 
face and the full analytic expansion of a'^^{B) in the powers of 1/5 can be 
obtained for ljbt ^ I. As we will see later this situation does not take place 
even for more general cases of quasiperiodic topologically regular trajectories 
on the Fermi surface (as well as for more complicated "chaotic" trajecto- 
ries). However, the expressions for ij'^^{B) for the case of the open periodic 
trajectories give also the correct main approximation to a^''{B) in the main 
order of l/S as we will see below. The corresponding dependence was also 
assumed in the papers where the more complicated cases of open trajec- 
tories appeared. They were also used in [201 1211 IHI] where the " Topological 
quantum numbers" corresponding to general topologically regular trajecto- 
ries were introduced. As we will see, however, the "ergodic dynamics" of (jH) 
in the more general cases can be observed in the next orders in 1/i? in precise 
conductivity measurements. Let us give here the corresponding expressions 
for conductivity for cases of compact and open periodic trajectories obtained 
in pp. We assume here that the z-axis is always directed along the magnetic 
field B and the x-axis in the plane n(B) (orthogonal to B) is directed along 
the mean direction of the periodic trajectory considered in p-space. The 
asymptotic values of cr*^ can then be written in the following for in the main 
order on 1/B ([I4): 

Case 1 (compact orbits): 



Here ~ means "of the same order in 00 bt and * are some constants ~ 1. 




(14) 



Case 2 (open periodic orbits): 




(15) 
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Let us mention also that the relations (fT^ - (fT3j) give only the absolute values 
of a'^. 

Let us mention also that the electron dynamics in x-space can be de- 
scribed by additional system 

and can be easily reconstructed for any known trajectory of (@]) in p-space. 
Such the a;?/-projection of any trajectory in x-space can be obtained just by 
rotation of corresponding trajectory of (0} by n/2. Easy to see then that the 
mean direction of open periodic electron trajectory in x-space coincides with 
the ?/-direction in our coordinate system. The last fact can be clearly seen 
in the form (|15|) of the tensor o"*'^ as the strong anisotropy of longitudinal 
conductivity in the plane n(B) for ljbt ^ 1. 

Let us now consider the cases of Topological rank 2 and 3 of the Fermi 
surface. 

3) The Fermi surface has Topological Rank 2. 

It can be easily seen that this case gives much more possibilities for the 
existence of open orbits for different directions of the magnetic field. In 
particular, this is the first case where the open orbits can exist for the generic 
direction of B with irrationality 3. So, in this case we can have the whole 
regions on the unit sphere such that the open orbits present for any direction 
of B belonging to the corresponding region. It is easy to see, however, that 
the open orbits have also a quite simple description in this case. Namely, 
any open orbit (if they exist) lies in this case in the straight strip of the finite 
width for any direction of B not orthogonal to the integral planes given by 
the components of Rank 2. The boundaries of the corresponding strips in the 
planes n(B) orthogonal to B will be given by the intersection of n(B) with 
the pairs of integral planes bounding the corresponding components of Rank 
2. It can be also shown (^7], P^) that every open orbit passes through the 
strip from —oo to -|-oo and can not turn back. 

According to the remarks above the contribution to the conductivity given 
by every family of orbits with the same mean direction reveals the strong 
anisotropy when u bt —>■ oo and coincides in the main order with the formula 
(fT^ for the open periodic trajectories. 

For purely rational directions of B we can have the situation when the 
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open trajectories with different mean directions present on different compo- 
nents of the Fermi surface. For example, for the "exotic" surface shown at 
Fig. IHl b we will have the periodic trajectories along both the x and y di- 
rections in different planes orthogonal to B if B is directed along the z axis. 
However, it can be shown that for any direction of B which is not purely 
rational this situation is impossible. We have so that for any direction of B 
with irrationality 2 or 3 all the open orbits will have the same mean direction 
and can exist only on the components of Rank 2 with the same (parallel) in- 
tegral orientation. This statement is a corollary of more general topological 
theorem which we will discuss below. 

At last we note that the directions of B orthogonal to one of the compo- 
nents of Rank 2 are purely rational and all the non-singular open orbits (if 
they exist) are rational periodic in this case. For any family of such orbits 
with the same mean direction the corresponding contribution to the conduc- 
tivity can then be written in the form (fT3jl in the appropriate coordinate 
system. However, the direction of open orbits can not be predicted apriori 
in this case. 

Let us say that the trajectories of this type have already all the features 
of the general topologically integrable situation and the topologically regular 
trajectories appearing in the more complicated case of topological rank 3 
have in fact very similar structure to the described above. However, we 
will need the special procedure of topological reconstruction preserving all 
the open trajectories on the Fermi surface which permits to see this fact in 
more complicated general situation. We will consider here in details also the 
"ergodic" properties of such trajectories on the corresponding pieces of the 
Fermi surface using the topological structure of this kind. 

Let us start now with the most general and complicated case of arbitrary 
Fermi surface of Topological rank 3. 

We describe first the convenient procedure ([221, EI] ) reconstruction of 
the constant energy surface when the direction of B is fixed. 

As we said already, we admit only the non-degenerate singularities of 
the dynamical system (@)) having the form of the non-degenerate poles or 
non-degenerate saddle points. The singular trajectories passing through the 
critical points (and the critical points themselves) divide the set of trajecto- 
ries into the different parts corresponding to different types of trajectories on 
the Fermi surface. As we also already said, we will not be interested here in 
the geometry of compact electron orbit in the "geometric limit" ubt —>■ oo. 
It's not difficult to show that the pieces of the Fermi surface carrying the 
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Figure 9: The cylinder of compact trajectories bounded by the singular or- 
bits. (The simplest case of just one critical point on the singular trajectory.) 

compact orbits can be either infinite or finite cylinders in bounded by 
the singular trajectories (some of them maybe just points of minimum or 
maximum) at the bottom and at the top (see Fig. 

Easy to see that the first case corresponds then to the whole component 
of rank 1 carrying just the compact trajectories while the second case gives 
the part of more complicated Fermi surface filled by the compact trajectories. 

Let us remove now all the parts containing the non-singular compact 
trajectories from the Fermi surface. The remaining part 

S'i?/ (Compact Nonsingular Trajectories) = Uj S'j 

is a union of the 2-manifolds Sj with boundaries dSj who are the compact 
singular trajectories. The generic type in this case is a separatrix orbit with 
just one critical point like on the Fig. El 

Easy to see that the open orbit will not be affected at all by the construc- 
tion described above and the rest of the Fermi surface gives the same open 
orbits as all possible intersections with different planes orthogonal to B. 
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Definition 6. 

We call every piece Sj the "Carrier of open trajectories". The tra- 
jectory is "chaotic" if the genus g{Sj) is greater than 1. The case g{Sj) = 1 
we call " Topologically Completely Integrable" . ^ 

Let us fill in the holes by topological 2D discs lying in the planes orthog- 
onal to B and get the closed surfaces 

Sj = Sj U (2 — discs) 

(see Fig. 

This procedure gives again the periodic surface S^ after the reconstruction 
and we can define the "compactified carriers of open trajectories" both in 
and T^. 

Let us formulate now the main topological theorems concerning the ge- 
ometry of open trajectories which made a breakthrough in the theory of such 
dynamical systems on the Fermi surfaces (|16], [T^). 

Theorem 1. [TB] 

Let us fix the energy level S'e and any rational direction Bq such that no 
two saddle points on S^ are connected in by the singular electron trajectory. 
Then for all the directions of B close enough to Bq every open trajectory 
lies in the strip of the finite width between two parallel lines in the plane 
orthogonal to B. 

In fact, the proof of the Theorem 1 was based on the statement that genus 
of every compactified carrier of open orbits Sj is equal to 1 in this case. 

Theorem 2. ^ 

Let a generic dispersion relation 

e(p) : ^ R 

be given such that for level e(p) = eo the genus g of some carrier of open 
trajectories Si is greater than 1. Then there exists an open interval (61,62) 

■^Such systems on were discussed for example in ^S]; the generic open orbits are 
topologically equivalent to the straight lines. Ergodic properties of such systems indeed 
can be nontrivial as it was found by Ya. Sinai and K.Khanin in |^. 
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Figure 10: The reconstructed constant energy surface with removed compact 
orbits the two-dimensional discs attached to the singular orbits in the generic 
case of just one critical point on every singular orbit. 
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containing cq such that for all e ^ tQ in this interval the genus of carrier of 
open trajectories is less than g. 

Let us note now that the genus 1 corresponds precisely to the two- 
dimensional tori embedded to the three dimensional torus T'^ as the com- 
pactified carriers of open trajectories on the Fermi surface. The correspond- 
ing components in the covering space can thus be either the periodically 
deformed cylinders or the periodically deformed integral plane in filled by 
the open orbits. Easy to see then that the first case can arise only for the 
directions of B orthogonal to some integer vector in (i.e. to the mean 
direction of the corresponding component carrying open orbits) and corre- 
sponds to the periodic open orbits in n(B). The second case corresponds 
to the more general situation and can arise for generic directions of B being 
stable with respect to all small rotations of the magnetic field. Easy to see 
the direct analogy of this situation with the cases of Fermi surfaces of rank 1 
and 2 respectively. Actually, we show in this construction that for the fixed 
direction of B we can replace the Fermi surface by some other surface having 
the simpler structure which gives us the same open trajectories in all planes 
n(B) orthogonal to the magnetic field. As can be easily seen also, all the 
trajectories can correspond only to topologically integrable case if the carri- 
ers of open trajectories all have the genus 1. The mean direction of generic 
open trajectories can be obtained then as the intersections of the planes or- 
thogonal to B with the integral planes in reciprocal lattice represented by 
the periodically deformed integral planes carrying open orbits. 

The Theorem 2 claims then that only the " Topologically Integrable case" 
can be stable with respect to the small variations of energy level and has a 
generic properties in this situation. It follows also from Theorem 1 that this 
case is the only stable case with respect to the small rotations of B and then 
should be considered as the generic situation among all the situations when 
the non-compact electron trajectories arise on the Fermi surface. 

The important property of the compactified components of genus 1 arising 
for the generic directions of B is following: they are all parallel in average in 
M.^ and do not intersect each other. This property, playing the crucial role 
for conductivity phenomena, was first pointed out in [20.J and called later the 
"Topological resonance" for the system According to this property, all 
the stable topologically regular open orbits in all planes orthogonal to B have 
the same mean direction and then give the same form ()15|) of contribution to 
conductivity in the appropriate coordinate system common for all of them. 
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This fact gives the experimental possibihty to measure the mean direction 
of non-compact topologically regular orbits both in x and p spaces from the 
anisotropy of conductivity tensor a*'^. From the physical point of view, all 
the regions on the unit sphere (giving the directions of magnetic field) where 
the stable open orbits exist can be represented as the "stability zones" Qq 
such that each zone corresponds to some integral plane Fq, common to all 
the points of stability zone Q^- The plane Fq, is then the integral plane in 
reciprocal lattice given by the integral mean direction of the components of 
genus 1 carrying open orbits and defining the mean directions of open orbits 
in p-space for any direction of B belonging to Qa just as the intersection 
with the plane orthogonal to B. As can be easily seen from the form of (|15|) 
this direction always coincides with the unique direction in corresponding 
to the decreasing of longitudinal conductivity as — >• cxo. 

The corresponding integral planes Fq, can then be given by three integer 
numbers {nl^jn^jnl^) (up to the common multiplier) from the equation 

nl[x\i + nllK]2 + nl[x\3 = 

where [x]j are the coordinates in the basis {gi, g2, gs} of the reciprocal lattice, 
or equivalently 

n^(x, li) + (x, I2) + nli^, I3) = 

where {li, I2, 13} is the basis of the initial lattice in the coordinate space. 

We see then that the direction of conductivity decreasing 17 = {rji,rj2,rj^) 
satisfies to relation 

nlifi, li)+n2(77, h)+nl{fi,h) = Q 

for all the points of stability zone VL^ which makes possible the experimental 
observation of numbers (ra^, ra^, n^). 

The numbers nl^-, n^a) were called in [201 the "Topological Quantum 
numbers" of a dispersion relation in metal. 

Let us add also that the number of tori T? being even can still be dif- 
ferent for the different points of stability zone ^q,. We can then introduce 
in the general situation the "sub-boundaries" of the stability zone which are 
the piecewise smooth curves inside Vl^ where the number of tori generically 
changes by 2. The asymptotic behavior of conductivity will still be described 
by the formula (fTK|) in this case but the dimensionless coefficients will then 
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"jump" on the sub-boundaries of stability zone. Let us however mention 
here that this situation can be observed only for rather complicated Fermi 
surfaces. 

Let us say now some words about the special situations when the carriers 
of open trajectories have genus more than 1 which can also happen for rather 
complicated Fermi surfaces. 

It was first shown by S.P.Tsarev ([2H]) the more complicated chaotic 
open orbits can still exist on rather complicated Fermi surfaces Sp- Such, the 
example of open trajectory which does not lie in any finite strip of finite width 
was constructed. The corresponding direction of B had the irrationality 2 
in this example and the closure of the open orbit was a "half of the surface 
of genus 3 separated by the singular closed trajectory non-homotopic to zero 
in T'^. However, the trajectory had in this case the asymptotic direction 
even not being restricted by any straight strip of finite width in the plane 
orthogonal to B. 

As was shown later in [221 this situation always takes place for any chaotic 
trajectory for the directions of B with irrationality 2. We have so, that for 
non-generic "partly rational" directions of B the chaotic behavior is still 
not "very complicated" and resembles some features of stable open electron 
trajectories. 

The corresponding asymptotic behavior of conductivity should reveal also 
the strong anisotropy properties in the plane orthogonal to B although the 
exact form of a^^ will be slightly different from (jl5|) for this type of trajecto- 
ries. By the same reason, the asymptotic direction of orbit can be measured 
experimentally in this the direction of lowest longitudinal conduc- 

tivity in according to kinetic theory. The measure of the corresponding 
set on the unit sphere is obviously zero for such type of trajectories being 
restricted by the measure of directions of irrationality 2. We will consider 
here the Tsarev example in more details below and discuss the corresponding 
conductivity behavior for B ^ oo. 

The more complicated examples of chaotic open orbits were constructed in 
[221 for the Fermi surface having genus 3. The direction of the magnetic field 
has the irrationality 3 in this case and the closure of the chaotic trajectory 
covers the whole Fermi surface in T^. These types of the open orbits do not 
have any asymptotic direction in the planes orthogonal to B and have rather 
complicated form "walking everywhere" in these planes. Let us also discuss 
later this case in more details. 
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Let us formulate now the recent topological results concerning the general 
behavior of quasiclassical electron orbits on the energy levels of the generic 
dispersion relation e = e(p) on T'^. 

The systematic investigation of the open orbits was completed after the 
works Uni US I2ni by I.A.Dynnikov (see In particular the total 

picture of different types of the open orbits for generic dispersion relations 
was presented in j^T]. Let us just formulate here the main results of 177] 
in the form of Theorem. 

Theorem 3 ([22], [2I|)- 

Let us fix the dispersion relation e = e(p) and the direction o/B of irra- 
tionality 3 and consider all the energy levels for e^m < e < ^max- Then: 

1 ) The open electron trajectories exist for all the energy values e belong- 
ing to the closed connected energy interval ei(B) < e < e2(B) which can 
degenerate to just one energy level ei(B) = e2(B) = eo(B). 

2) For the case of the nontrivial energy interval the set of compactified 
carriers of open trajectories 5*^ is always a disjoint union of two-dimensional 
tori in for all ei (B) < e < €2 (B) . All the tori for all the energy levels 
do not intersect each other and have the same (up to the sign) indivisible 
homology class c G H2{T^,Z), c 7^ 0. The number of tori is even for 
every fixed energy level and the corresponding covering in is a locally 
stable family of parallel ("warped") integral planes C with common 
direction given by c. The form of described above is locally stable with 
the same homology class c G 7^2 (T^) under small rotations of B. All the 
open electron trajectories at all the energy levels lie in the strips of finite 
width with the same direction and pass through them. The mean direction of 
the trajectories is given by the intersections of planes n(B) with the integral 
family for the corresponding "stability zone" on the unit sphere. 

3) The functions ei(B), e2(B) defined for the directions of B of irra- 
tionality 3 can be continuated on the unit sphere 5*^ as the piecewise smooth 
functions such that ei(B) > e2(B) everywhere on the unit sphere. 

4 ) For the case of trivial energy interval ci = 62 = eo the corresponding 
open trajectories may be chaotic. Carrier of the chaotic open trajectory is 
homologous to zero in H2{T^,1j) and has genus > 3. For the generic energy 
level e = Co the corresponding directions of magnetic fields belong to the 
countable union of the codimension 1 subsets. Therefore a measure of this 
set is equal to zero on S"^. 
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The whole manifold is always homologous to zero in and all the 
two dimensional tori can be always divided in two equal groups {T^_^}, 
{¥?_} according to the direction of the electron motion. As can be proved 
using Theorem 1 (|2Z|) the "stability zones" form the everywhere dense set 
on the unit sphere for the generic dispersion relations. All the non-compact 
trajectories stable under the small rotations of B should have thus the form 
described above. 

Let us however remind here that for the conductivity phenomena in metals 
only the Fermi level e(p) = ep is important and we can not get any informa- 
tion about stability zones of dispersion relation such that ep ^ [^2(6), ei(B)]. 

Let us say also that Theorem 3 describes the generic situation of the 
directions of B of irrationality 3 and some additional features can arise for 
the purely (irrationality 1) or "partly" (irrationality 2) rational directions 
of B. We will not discuss here all these features in details and just make a 
reference on the articles Ell where the general situation was considered. 

4 The behavior of conductivity in the case of 
topologically regular open orbits. 

Let us now discuss briefly the conductivity behavior in "geometric strong 
magnetic field limit" . We prove below the analog of Lifshitz-Azbel-Kaganov- 
Peschanskii formulae in this case ([HIS IS])- We will start with generic "topo- 
logically regular" open orbits having the "strong asymptotic directions in M^. 

Let us come back to the formula (fT^ and consider the case of the open 
orbits having the form shown at Fig. El a. Let us take the x-axis in 
with the direction corresponding to the mean directions of open orbits in 
p-space and the z-axis directed along the magnetic field B. The projection 
of trajectory on the plane orthogonal to B will have the same form in x- 
space just rotated by 7r/2 such that the mean direction of open orbits will be 
directed along the y-axis in our coordinate system. 

We can write then that generally 

K) = , K) ^ , i^vi;) ^ 

on these trajectories, where (. . .) means the averaging along the trajectory 
according to system (|H). 

Moreover we can claim that the function 
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eB /■* 

p^(p,-t)-p^(p,o) =— vi.{jp,-t')dt' 

is a bounded function for any trajectory shown at Fig. IHl a, as follows from 
its geometric form. 

Let us consider now the components a*^ (or cr^*) of the conductivity tensor 
(j**^(B). Using the integration by parts we can write: 



eBr 



dfr 
de 



V(P)x 



X 



r+oo 

/ {p'{p,-t)-p\p,0))e''/^dt 
Jo 



{2Tihf 
(16) 



For the component a we can use then the formula 



a-(B) 



{eBfr 



T 



de 



i dipyip,~t)-pyip,o)f 

2 dt 



-"^dt 



d^p 

{27C hy 



.eBr 
The integral 



d^p 



r+oo T 

1 = ^{p'{p,-t)-p'ip,0)fe-'/^dt 

can be evaluated as J < 1/2tPq where Pq is the common constant bounding 
the values |p^(p, —t) — p^(p, 0)| for all the regular trajectories (the width of 
straight strip). The value ct^^(B) can then be written as 



a^\B) - + o 



(17) 



when B — >■ oo. 
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The value of constant a can then be defined in terms of the integral 



over the Fermi surface and is proportional to the square of effective width of 
straight lines bounding the topologically regular trajectories. 

Let us say now that the total expansion in (1/-B) can not be actually 
made for the generic regular trajectories because of their ergodic behavior. 
However, the difference between the quasiperiodic and periodic trajectories 
does not appear in the main terms of the formula (fT3j) and can be detected 
only in the next smaller approximations for B ^ oo. This fact permits 
actually to use the formula (fT^ also for regular quasiperiodic trajectories 
when speaking just about the geometric properties of a^^ . Let us add also 
that the next approximations to (fT3jl depend on whether the saddle critical 
points really present on the carriers of open orbits or not. 

For the components o"*^ or o"^*, i = 2, 3 we can not use the second inte- 
gration by parts and we should put then 



according to the formula (fTB|) . 

Again the same remarks about next approximations in o"*^ and cr^* can 
also be made in this case. 

Let us consider now components a^^ where both i,j 7^ 1. We can write 
the the formula (fT^ using the averaged values of f and Vg^ in the following 
way: 





(18) 






(19) 
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(It's not difficult to see the the "cross terms" will disappear after the inte- 
gration over T^.) 

The ffist part of this formula gives the finite values for cr*^(B) when 
B ^ oo while the second becomes zero in the same limit. Easy to see that 
the last statement is just a simple corollary of the fact that the integral with 
decreasing exponent approaches the mean value on the trajectory in the limit 
B ^ oo (up to the factor 1/r) according to the system (jH). Let us say again, 
however, that the -B-dependence is more complicated in general here than in 
the case of purely periodic trajectories and can not be generically expanded 
in the powers of 1/B. 

Thus we can see that the formula (|15p gives the main part of (T*'^(B) also 
for general topologically regular trajectories with irrational mean directions. 
The cases of compact and purely periodic trajectories, however, are charac- 
terized by the existence of full analytic expansions in 1/B which do exist in 
general for more complicated quasiperiodic orbits. 

Let us make also some additional remark on the periodic and quasiperi- 
odic open trajectories. Namely, the periodic open orbits arise actually in 
every "stability zone" when the intersection of the carrier of open orbits 
with the plane n(B) has the rational mean direction. Easy to see that this 
situation can appear if the direction of B has irrationality 2 and always arise 
for the purely rational directions of B. This actually means that the corre- 
sponding orbits in are not everywhere dense on the tori (defined for 
a given stability zone) but become closed (in T^) unlike the situation with 
irrational mean direction of open orbits. The mean values of v^^. then ac- 
tually depend on the open periodic trajectory and do not coincide with the 
value of Vg^ averaged over the carrier of open trajectory. We know just that 
averaging of these mean values over all the periodic trajectories will give us 
the values of (Vgr) close to the same values for close irrational directions of 
B. However, the limit B ^ oo for tensor ct*^(B) {i,k = 2,3) will be given 
by formula 



containing the integration of products (f*r)(p)(fgr)(p) which can not be re- 
placed by the "global" mean values (Vgr), {Vgr)- As a corollary the corre- 
sponding values of a^^ will be actually different in this limit from the case 
of purely irrational directions of B. Such we can claim that the longitudinal 
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conductivities cr^^ and cr^^ should be actually bigger for B — oo than the 
same values for close irrational directions of B. 

We can see then that except the smooth dependence on the direction 
of B (connected with the dependence on B of total phase volume of open 
trajectories) we will have also the "sharp" peaks in the conductivity for the 
rational mean directions of open orbits. The corresponding "peaks" however 
will be quite small for the rational directions with big denominators. Easy to 
see also that this behavior does not affect the "Topological characteristics" 
of conductivity introduced above. 

5 The chaotic cases. 

Let us say now some words about chaotic trajectories which can arise in the 
special cases for rather complicated Fermi surfaces. We will first describe 
the Tsarev's example of chaotic trajectory having an asymptotic direction 
in M.^ (EHl)- Let us consider the Fermi surface Sp consisting of two sets of 
horizontal parallel layers {Q+} and {Q-} connected by two sets of inclined 
cylinders {Ci}, {C2} in (gee Fig. 

We will assume that both types of the cylinders have the identical forms 
and two cylinders shown at Fig. E]have the same vertical symmetry plane 
n which intersects the horizontal layers in some irrational direction a (Fig. 

We can assume now that the dark region on Fig. IT^ represents actually 
the "jumping place" where the intersection curves 11' fl Sp "jump" from one 
horizontal layer to another (Fig. [TT|) for any plane 11' parallel to 11. Let us 
put now n = n(B) and consider the intersections 11 fl S'i;' as the electron 
trajectories in p-space for given B. 

We will put for simplicity that the horizontal periods of Fermi surface are 
orthogonal to each other and we can divide the horizontal layers into identical 
periodic strips parallel to one of the periods (say vertical strips formed by 
the repeated domains shown at Fig. IT^ . Easy to see that for our specific 
direction a any trajectory will either pass through the strip being at the same 
layer in or meet the pair of cylinders (of identical form) and jump by two 
layers preserving the layer typer ("positive" or "negative"). The horizontal 
projection of corresponding orbits parts are always the same in this situation 
while the vertical displacements can be either zero or the vertical period of Sp 
if the vertical jump takes place in the corresponding strip. The probability 
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Figure 11: The Fermi surface for Tsarev's example of chaotic trajectory. 
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Figure 12: The intersection of the horizontal planes by the vertical plane with 
irrational direction passing through the axes of inclined cylinders in Tsarev 
example. 
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of jump is proportional to the ratio 6/ a (Fig. IT^ and all the non-singular 
trajectories have then the same asymptotic directions. 

It can be also easily seen that the trajectories belonging to different types 
of layers go and "jump" in the opposite directions being divided by the 
singular trajectories passing through the saddle-points on the sides of the 
cylinders Ci and C2. 

It can be proved, however, that there exist the special values of a such 
that the trajectories will not belong here to the straight strips of any finite 
width as follows from the numbers theory. In this situation we can not use 
anymore the same arguments which we used for calculating of cr^* and cr^* in 
the case of topologically regular trajectories and so we don't have formulas 
(II 7|) and (jl8|) in this situation. However, we can still use the fact 



(in the appropriate coordinate system) also in this case. Using the same 
arguments for the formula (|T^ we can prove then the following formulae for 
the behavior of (T*^(B): 



which replaces the formula (|T3j) for the case of Tsarev's ergodic trajectories. 
Let us omit here all the details of Tsarev's chaotic trajectories and just point 
out that the asymptotic direction of ergodic trajectory of Tsarev type can 
be also observed experimentally due to the same reasons as in the case of 
topologically regular trajectories. However, unlike the topologically regular 
case, the ergodic trajectories of Tsarev type are unstable with respect to 
generic small rotations of B and do not correspond to "stability zones" on 
the unit sphere. 

Let us say now some words about more general ergodic trajectories of 
Dynnikov type which do not have any asymptotic direction in M^. We will 
not describe here the corresponding construction (see j^) and just give the 
main features of such trajectories. 

First of all, these trajectories can arise only in the case of magnetic field 
of irrationality 3 and the corresponding carriers have then the genus > 3. 
This kind of trajectories are completely unstable with respect to the small 
rotations of B and can be observed just for special fixed directions of B 



Kr) = , {vy,) ^ , ivi;) ^ 




(20) 
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in the case of rather comphcated Fermi surfaces. The approximate form of 
trajectories of this kind is shown at Fig. |2l b. The carrier of such trajectory 
can then resemble the Fermi surface shown at the Fig. d and is homologous 
to zero in the 3-dimensional torus T^. Moreover, if the genus of the Fermi 
surface is not very high (< 6) it can always be stated that the corresponding 
carrier of open ergodic trajectories is invariant under the involution p — p 
(after the appropriate choice of the initial point in T^). The ergodicity of the 
open trajectories on the carrier gives then immediately the relations: 



for all three components of the group velocity on any of such trajectories. 
This important fact leads to the rather non-trivial behavior of corresponding 
contribution to the conductivity for B ^ oo. Namely, using the same formula 
fll2|) in this case we can show that all the components of the corresponding 
contribution to (T*'^(B) become actually zero in the limit B ^ oo (|24j). We 
can write then for this contribution: 



for B oo.^ 

We see then that the chaotic trajectories of Dynnikov type do not give any 
contribution even for conductivity along the magnetic field B for rather big 
values of B. In the work [25 also the special "scaling" asymptotic behavior 
of cr*'^(B) were suggested. Let us note, however, that the full conductivity 
tensor include also the contribution of compact (closed) electron trajectories 
having the form ()14p which presents in general as the additional contribution 
in all the cases described above. We can so claim that the ergodic behavior of 
Dynnikov type does not actually completely removes the conductivity along 
the magnetic field B because of the contribution of compact trajectories. 
However, the sharp local minimum in this conductivity can still be observed 
in this case since a part of the Fermi surface will be effectively excluded from 
the conductivity in this situation. 

^Actually the component cr^^(B) contains the non- vanishing term of order of T^/e|, 
for B ^ oo for non-zero temperatures ( [21) ) . However, this parameter is very smah for 
the normal metals and we don't take it here in the account. 



= , «) = , «) = 




(21) 
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It can be proved (see j2Z|) that for generic Fermi surfaces the measure of 
directions of magnetic field B where the chaotic behavior of Dynnikov type 
can be found on the Fermi surface is zero. Let us say that the restriction 
on just one energy level is connected closely with the situation in metals 
where only the energy levels close to are important. Let us formulate also 
the more general conjecture of S.P.Novikov about the total set of "chaotic 
directions" on the unit sphere, where we do not restrict the consideration to 
just one energy level: 

Conjecture. (S.P.Novikov). 

The total set of the directions of B corresponding to the chaotic behavior 
has the measure for the whole generic dispersion relation and the Hausdorf 
dimension strictly less than 2. 

Let us mention also that in the paper the possibilities of the investiga- 
tion of total topological characteristics of whole dispersion relation e(p) were 
discussed. In particular the behavior of electrons injected in the empty band 
of semiconductor in the presence of magnetic field was considered. However, 
the corresponding magnetic fields should be extremely high (~ lO^T/) in 
this case which make such experiments very difficult. Most probably this 
situation should be considered now just as theoretical possibility. 

6 The different regimes of conductivity be- 
havior. Classification. 

We are going to give now the full classification of possible conductivity 
regimes corresponding to different topological types of trajectories given by 
system (0} . First we will need to add here the non-generic cases of irrational- 
ity 1 and 2. The Theorem 3 should be slightly modified in this case but has 
the same main features as in the case of fully irrational magnetic field f|27j). 
Namely, the set of carriers of open trajectories can contain now the two- 
dimensional tori T^' having the zero homology class in H2{T^) in addition to 
the family of parallel tori with non-zero homology classes described above. 
The corresponding covering of these components in are "warped" peri- 
odic cylinders and all the open trajectories belonging to these components 
are purely periodic. As it is easy to see these components of are stable 
with respect to the small rotations of B in the plane orthogonal to the axis 
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of cylinder and disappear after any other small rotation. The part (1) of 
the Theorem 2 will be true also for rational or "partly rational" directions 
of B with some connected energy interval e'^(B) < e < e2(B). However, the 
boundary values e'^(B), e2(B) do not necessarily coincide in this case with 
the values of piecewise smooth functions ei(B), e2(B) defined everywhere on 
5*^ according to Theorem 3 ( [2Z| ) ■ Namely, we will have instead the relations 
— ^i(-^) — ^2(B) < e2(-B) such directions of B where all the 

components of belonging to intervals [e'^(B), ei(B)), (e2(B), e2(B)] consist 
of the tori homologous to zero in T^. As we will see this cases can be ob- 
served experimentally for those B where the Fermi level lies in the one of 
such intervals and only the "partly stable" non-compact trajectories exist on 
the Fermi surface. 

The "partly-stable" cylinders described above do not intersect the "ab- 
solutely stable" components of and all the open trajectories will still have 
the same mean direction if B is not orthogonal to corresponding integral 
plane Fq, (let us mention that all the trajectories lying on "regular" paral- 
lel integral planes in 5*^ C will be also periodic with the same period in 
this case) . The form of the conductivity tensor will still be described by the 
formula (fT3j) but the numerical values of dimensionless coefficients will jump 
for those non-generic directions of B where the situation described appears. 
In particular, the asymptotic behavior (fl^j) will arise on the one-dimensional 
curves for those directions of B where ep G [e[(B) , eiCB)) U (e2(B), e2(B)] 
being stable only for rotations of B in the corresponding direction. As fol- 
lows from the statements above the corresponding directions of B can have 
at most the irrationality 2 and the corresponding one-dimensional curves 
are always the parts of the circles orthogonal to some integer vector in the 
reciprocal lattice F*. 

Let us make now a special remark about the "Special directions" of B 
orthogonal to the integral planes Fq, if this direction belongs to the corre- 
sponding stability zone Qa- The direction of B is then purely rational and all 
the corresponding open orbits (if they exist) should be periodic in R^. How- 
ever, the mean directions of these open orbits can be different in this case for 
the different planes orthogonal to the magnetic field and the corresponding 
contributions to conductivity can not be then written in the form (fT3jl in the 
same coordinate system. The conductivity tensor a*^ can have then the full 
rank in the limit Br oo and the conductivity remains constant for all the 
directions in R^ in the strong magnetic field limit. This situation, however, 
is completely unstable and disappear after any small rotation of B. 
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Figure 13: The "singular periodic net" of open trajectories in the plane 
orthogonal to B. 

The second possibility in the case of such directions is that all the open 
trajectories become singular and form the "singular periodic nets" in the 
planes orthogonal to B (see Fig. IT^ . The asymptotic behavior of conduc- 
tivity is described then by formula (fT^ but is also completely unstable and 
changes to ()15|) after any small rotation of B. 

Let us add that the same situations can arise also on the stable two- 
dimensional tori Tf in this case where the directions of open orbits can not 
be defined anymore as the intersection of with the plane n(B). Such, we 
can have either the "singular nets" or the regular periodic open orbits on 
these tori for this special direction. Also we can have the open orbits with 
different mean directions on different tori in this case but the number of such 
tori should be > 4 for any physical type of dispersion relations. This situation 
can thus also be observed only for rather complicated Fermi surfaces. 
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Figure 14: The wide finite strip in the plane n(B) containing the open orbits 
for B close to "Special rational direction" within the stability zone. 

Let us mention also that for the directions of B close to this special 
rational directions the widths of the straight strips containing the regular 
open orbit can become very big in the planes n(B) (see Fig. IT^ . 

So from physical point of view the conductivity phenomena will not "feel" 
the mean directions of open orbits for the directions of B close enough to these 
special points even for rather big (but finite) values of magnetic field. Instead, 
the oscillations of the trajectory within the strip (Fig. IT^ will be essential for 
conductivity up to the rather big values of B such that ubt ~ L/pq (where 
L is the width of the strip and Pq is the size of the Brillouin zone). These 
situation exists, however, if the open orbits with corresponding direction exist 
also for B = Bq where Bq is the special rational direction in the stability 
zone. For B close to Bq we will observe then exactly this direction up to 
the values of B such that ubt ~ L/pq and then the regime will change to 
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the common situation corresponding to given stability zone. Experimentally 
we will observe then the "small spots" around these directions on the unit 
sphere where the anisotropy of a^^ corresponds to the direction of orbits for 
B = Bq. In the most complicated case when we have the open orbits with 
different mean directions for B = Bq we will have then the finite conductivity 
for all the directions in in the corresponding spot for rather big values of 
B. 

For the "special rational directions" corresponding to the "singular net" 
on the stable tori T^ the behavior of o"*'^ will correspond to the common 
form for a given stability zone even for B very close to Bq. However, the 
measure of open orbits will tend to zero as B — Bq in the stability zone. 
The dimensionless coefficients (*) in the formula (|15|) will vanish then for 
B — > Bq although the integral plane Fq, will be observable up to B = Bq. 

We mention here at last that both cases when Bq belongs or does not 
belong to the corresponding stability zone Qa are possible in the examples. 

We can describe now the total picture for the angle diagram of conduc- 
tivity in normal metal in the case of geometric strong magnetic field limit. 
Namely, we can observe the following objects on the unit sphere parameter- 
izing the directions of B: 

1) The "stability zones" Qa corresponding to some integral planes F^ in 
the reciprocal lattice ("Topological Quantum numbers"). This "Topological 
Type" of open trajectories is stable with respect to small rotations of B and 
this is the only open orbits regime which can have a non-zero measure on the 
unit sphere. All the "stability zones" have the piecewise smooth boundaries 
on S"^ and are given by the condition 

ei(B) < < e2(B) 

where ei(B), 62 (B) are the piecewise smooth functions defined in the Theo- 
rem 3. Generally speaking, this set is not everywhere dense anymore being 
just a subset of the corresponding set for the whole dispersion relation and 
have some rather complicated geometry on the unit sphere. 

The corresponding behavior of conductivity is described by the formula 
(lisp and reveals the strong anisotropy in the planes orthogonal to the mag- 
netic field. For rather complicated Fermi surfaces we can observe also the 
"sub-boundaries" of the stability zones where the coefficients in (fT3jl have 
the sharp "jump". 
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2) The net of the one- dimensional curves containing directions of irra- 
tionahty < 2 where the additional two-dimensional tori (homologous to zero 
in T^) can appear. The corresponding parts of the net are always the parts 
of the big (passing through the center of S"^) circles orthogonal to some re- 
ciprocal lattice vector where the condition 

e[{B) < ep < e',{B) 

is satisfied. The asymptotic behavior of conductivity is given again by the 
formula (fT3j) . 

Let note also that these special curves on S*^ can be considered actually 
as the reminiscent of the bigger stability zones if we don't restrict ourselves 
just by one Fermi surface. The structure of such sets thus can be used to 
get more information about the corresponding total structure for the whole 
dispersion relation in metal. In particular, the mean direction of such open 
orbits coincides with the mean direction of the generic open orbits in the 
intersections of the net with the stability zones (except the " Special rational 
directions"). The corresponding conductivity tensor is given then by the 
same formula (|15p where the dimensionless coefficients "jump" on the curves 
of the net. 

3) The "Special rational directions". 

Let us remind that we call the special rational direction the direction of 
B orthogonal to some plane in case when this direction belongs to the 
same stability zone on the unit sphere. We can have here all the possibilities 
described earlier for this situation (i.e. regular behavior with vanishing coef- 
ficients in (fT3j) . spots with isotropic or anisotropic behavior of conductivity 
different from the given by corresponding "Topological quantum numbers", 
"partly stable" isotropic or anisotropic addition to the (jl5|) . etc.) 

4) The chaotic open orbits of Tsarev type (B of irrationality 2). 

We can have points on the unit sphere where the open orbits are chaotic 
in Tsarev sense. All open trajectories still have the asymptotic direction 
in this case and the conductivity reveals the strong anisotropy in the plane 
orthogonal to B as i? ^ oo. The B dependence, however is slightly different 
from the formula (jl5|) in this case. 

5) The chaotic open orbits of Dynnikov type (B of irrationality 3). 
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For some points on S'^ we can have the chaotic open orbits of Dynnikov 
type on the Fermi surface. At these points the local minimum of conductivity 
along the magnetic field is expected. The conductivity along B however 
remains finite as 5 — > cxd in general situation because of the contribution of 
compact trajectories. 

6) At last we can have the open regions on the unit sphere where only the 
compact trajectories on the Fermi level are present. The asymptotic behavior 
of conductivity tensor is given then by the formula ()14|) . 

Let now point out some new features connected with the "magnetic break- 
down" (self-intersecting Fermi surfaces) which can be observed for rather 
strong magnetic fields. Up to this point it has been assumed throughout that 
different parts of the Fermi surface do not intersect with each other. How- 
ever, it is possible for some special lattices that the different components of 
the Fermi surface (parts corresponding to different conductivity bands) come 
very close to each other and may have an effective "reconstruction" as a result 
of the "magnetic breakdown" in strong magnetic field limit. In this case we 
can have the situation of the electron motion on the self-intersecting Fermi 
surface such that the intersections with other pieces do not affect at all the 
motion on one component. (The physical conditions for the corresponding 
values of B can be found in [ID])- In this case the picture described above 
should be considered independently for all the non-selfintersecting pieces of 
Fermi surface and we can have simultaneously several independent angle dia- 
grams of this form on the unit sphere. Such we can have here the overlapping 
stability zones where the open orbits can have different mean directions. The 
correspondent conductivity tensor will then be given just as a sum of all con- 
ductivity tensors corresponding to different non-selfintersecting components. 
(The problem of the magnetic breakdown was brought to the authors' atten- 
tion by M.I.Kaganov.) 
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